Two dimensional dipolar scattering with a tilt 
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We study two body dipolar scattering in two dimensions with a tilted polarization axis. This tilt 
reintroduces the anisotropic interaction in a controllable manner. As a function of this polarization 
angle we present the scattering results in both the threshold and semi-classical regimes. We find a 
series of resonances as a function of the angle which allows the scattering to be tuned. However the 
character of the resonances varies strongly as a function angle. Additionally we study the properties 
of the molecular bound states as a function of the polarization angle. 

PACS numbers: 34.20.Cf,34.50.-s,34.20.-b 



I. INTRODUCTION 

There are exciting proposals based on dipolar gases 
in two-dimensional (2D) geometries, such theories show 
dipolar systems will lead to exotic and highly correlated 
quantum systems [lj. There has also been tremendous 
advances in the production of ultracold polar molecules 
Q, especially at JILA where a group has produced ul- 
tracold RbK in ID optical lattice Q. The group used 
an electric field to align the molecular dipole moments 
(d) along the z direction perpendicular to the plane of 
motion [p). The tight trapping geometry and the dipo- 
lar interaction were used to inhibit the molecules from 
reaching their short range interaction where they would 
chemically react 

There are alternative molecular systems which will not 
chemically react and less restrictive configurations can 
be considered. For example RbCs and NaCs are chemi- 
cally stable Q . One interesting possibility to control the 
properties of the these dipolar gases is to tilt the polar- 
ization axis into the plane of motion. Such a scenario 
has been considered by many-body theories: for example 
anisotropic superfluidity has been predicted @ , 2D dipo- 
lar fermions have been studied Q , and few body dipolar 
complexes have been investigated Q. However little is 
known about the nature of the scattering physics of such 
a 2D system with some in-plane polarization. That is the 
aim of this work; we study the scattering properties of the 
2D dipolar system when the polarization is not fixed out 
of the plane of motion. This reintroduces the anisotropy 
of the interaction controllably, such that there is a pre- 
ferred direction where the dipoles can line up in a head 
to tail fashion that is energetically favored to the side by 
side configuration, just as in 3D dipolar physics. But in 
this case the strength of the anisotropy can be controlled 
directly by the polarization angle. For small angles, there 
is only a weaker repulsion in one direction, but for large 
angles (near 7r/2) there is an attractive interaction. 

Some recent work has aimed at understanding the scat- 
tering behavior of dipoles in 2D and quasi- two dimensions 
(q2D). First, dipolar scattering 2D in the threshold and 
semi-classical regime was studied in pure 2D determin- 
ing the limiting behavior of such scattering @. Thenq2D 



was studied [10( and more recently Ref. [ll| introduced 
an elegant method to solve the full scattering problem. 
Other theories have focused on understanding scattering 
and chemical reactions in q2D [l^, EH and how to use 
the electric field and trap to control the scattering rate. 
There has also been some recent work on understand- 
ing the scattering and bound state structure of 2D layer 
dipolar systems [Til [l5| and in a layered system with a 
tilted polarization axis (l6| . 

In the next section we look at the basic scattering of 
the system and offer estimates of the scattering as a func- 
tion of the polarization angle. Then we look at the char- 
acter of the scattering resonances. Finally we study the 
molecules that can be formed, their binding energies, size, 
and shape as a function of polarization angle. 



II. EQUATIONS OF MOTION 

For this work we assume that the length scale of con- 
finement is much smaller than any other length. This 
effectively removes it from the scattering problem. Re- 
alistically this length scale will be important, but as a 
first study to provide useful estimates of the scattering 
this assumption is justified. The Schrodinger equation 
for two dipoles in 2D is 




Eijj. 



(1) 



where the dipoles are polarized along d = zcos(a) + 
x sin(a) with magnitude d and p = xx + yy- We solve 
this equation by using partial wave expansion: ip(p) = 
J2 m ( t>m(kp)e zm '^ I \fp. In this case the tilted polarization 
axis ruins the cylindrical symmetry, meaning that differ- 
ent azimuthal symmetries are coupled together. Impor- 
tant features of the interaction anisotropy can be distilled 
by looking at the matrix elements: 



(m|l - 3(d- pf\m!) = U„ 
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For a — 0, the system is totally repulsive and isotropic. 
Then as a is increased the isotropic repulsive term is 
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weakened in the x direction but it is still full strength 
in the y direction. This anisotropy enters as couplings 
between cha nnels with m and m ± 2. At the a c = 
sin _1 ('\/2/3) ~ 54.7 degrees or a/n ~ 0.3 there is no 
barrier to the short range and past this angle there is an 
attractive dipolar diagonal potential. 

Using the matrix elements in Eq. ([2]) and the dipo- 
lar length D = pd 2 /h 2 to rescale Eq. ([1}, we obtain a 
multi-channel radial Schrodinger equation describing 2D 
dipolar scattering with tilted polarization axis: 



Ydp 2 



1/4 



2p 2 
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where p — p/D and Ejj the dipolar energy is h 6 / p?d A . To 
perform the scattering calculation we start at po/D. We 
vary this parameter and use it to control the scattering 
properties of the system; po/D is used to tune the m=0 
scattering or the scattering length a/D. This distance 
signifies where the interaction becomes more complicated 
through transverse modes or system specific interactions 
becoming important. Thus a more sophisticated bound- 
ary condition is required at this wall. However for our 
initial treatment, we only demand the wave function be 
zero at po/D. 

It is worth while to comment that using Pq/D to pa- 
rameterize the scattering should be viewed as varying the 
electric field. D is pd 2 /h 2 and as and electric field is in- 
creased the induced dipole moment, d, becomes larger. 
Thus decreasing rhoo/D mimics an increasing electric 
field. Additionally, the correspondence of rhoo /D to a/D 
is unique but due to the nature of the problem (multi- 
channel) it is complex and numerically found. 

Before we present the full numerical scattering calcu- 
lations, we discuss the form of the free {JJ mm > = 0) 2D 
wavefunctions for both scattering and bound states. The 
scattering wavefunction is <p m (kp) = cos(S m ) f m (kp) — 
s'm(5 m )g m (kp) where S m is the scattering phase shift for 
the m partial wave and f m (g m ) is the regular (irregular) 
free solution. In 2D, it is y/kpj m (kp) (y/k~pN m ) where J m 
(N m ) is a Bessel (von Neumann) function and k — \J2pE. 
If the system were bound then the asymptotic wavefunc- 
tion for the m = is 4>b = y/KpKo(Kp) where Kq is the 
modified Bessel function and in the large p limit this de- 
cays as e~ Kp with k = y/—2[iEi, and Ej, is the binding 
energy. 

As in 3D, the scattering length is defined by when the 
zero energy wavefunction is zero, 'tjj(a) = </>o( a ) = 0: 
4>o(a) = cot(<5o)/o(a) — 9o( a ) = [13 ■ Then the scat- 
tering length can be computed with a = |ei cat(*o)-7 
where 7 is the Euler gamma function ~0.577. Conversely, 
the phase shift can be defined by the scattering length: 
cot(<5o) = -(ln(fca/2) + 7) as the first term in the effec- 
tive range expansion [l7|. This definition of the scatter- 
ing length is effectively energy independent once in the 
thresholds regime, Dk < 1. 




FIG. 1: (Color Online) The Born approximation (BA) for 
fcrmions shown as a function of a/n. We have plotted the 
energy independent: a/a'^ l jl(a = 0) as a function of a for 
Dk = 10" 1 (red open circle), 1CT 2 (blue +), 10~ 3 (green 
dash), 10 -4 (purple open diamond), and the BA is shown as 
a bold blue line. 



Using these wavefunctions to study the basic molecular 
properties, we start in the large a limit where the binding 
energy goes to zero. At moderate p, the wavefunction is 
essentially the same whether the two particles are a near 
zero energy scattering state or a loosely bound molecule. 
Using this fact, we match log-derivatives of the near-zero 
energy scattering wavefunction in terms of the scattering 
length from the short range to the long range asymptotic 
bound wavefunction, ^(ftp). Using the small argument 
expansions of both wavefunctions allows us to determine 
K in terms of a: n a = 2e~~ 1 /a and the binding energy 
follows: -Ah 2 e- 2 ^ /2pa 2 Using (j) b (n a p) will offer 

us many interesting analytic properties of the molecules 
in the large a limit. When compared to the full multi- 
channel numerical calculation, these analytic results pro- 
vide very good estimates of the molecular properties for 
large values of a/D. 



III. SCATTERING 

Now we look at the scattering properties of this system 
as a function of the polarization angle. To do this, we 
solve Eq. ([3]) with the Johnson Log-derivative propagator 
[l9| . We then extract the T- matrix, Tj f , which describes 
the scattering between channels i and /. The total cross 
section is a = -|X)i/l^/l 2 - The elastic cross section 
for m can be written as a m = ^sin 2 (5 m ), where 5 m is 
the scattering phaseshift for the m partial wave [2(| Hl[ . 
Sometimes, the most useful quantity is not the scatter- 
ing cross section, rather it is the dimensionless scattering 
rate ka. Plotting this quantity reveals the system in- 
dependent or universal scattering characteristics of the 
scattering dipolar system, as was shown in 3D by Refs. 
[22H24l | . We now present general trends of the 2D dipolar 
system with a tilted polarization axis. 
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In the threshold regime the Born approximation (BA) 
offers a good estimate of the scattering for non-zero par- 
tial waves 0, . This is most useful for estimating the 
cross section for identical fermions. In this model the 
dipoles are spinless and the way one models fermions or 
bosons is by imposing the symmetric or anti-symmetric 
requirement on the spatial wavefunction. This leads to 
the scattering properties of fermions being a sum of only 
the odd partial waves and bosons a sum of only the even. 
For the case of distinguishable particles, one sums up all 
of the partial waves. The BA result for this systems is: 



i 4(£>fc) 2 f 3 . . 2 

k °BA = (m 2_l)2 ^-2 Sln ^ 



(4) 
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There are two basic types of collisions given here: diago- 
nal and off-diagonal or m changing. For the off-diagonal 
scattering, there is a special case of p-wave collisions 
where ±1 goes to ^1 and has the same functional form 
as the diagonal contribution. Then for the other off- 
diagonal terms m m + 2, there is a distinct form. For 
identical fermions, the BA can be compared directly to 
the full scattering cross section with out worrying about a 
short range phase shift. This comparison is made in Fig. 
[T] where the full cross section is divided by the a = 
BA cross section. Plotting this removes the energy de- 
pendence of the cross section. The BA is shown as a 
thick blue line normalized by its a = value, and the 
full cross sections ares shown for Dk = 10~ 4 (violet open 
diamonds), 10~ 3 (dashed green), 10 -2 (black +), and 0.1 
(red open circles). Relating Dk to scattering energy is 
simply: E/E D = (Dk) 2 /2. 

In Fig. Q] the agreement is good for the full range 
of a, especially at small k. But it is worth noticing for 
Po/D = 0.01 there are two resonances as a/n goes to 1/2. 
This are most clearly seen in k = 0.1. For the smaller 
values of k they are narrow and the region where the as 
deviates from the BA are increasingly small. 

Alternatively, in the high energy regime we can esti- 
mate the cross section with the Eikonal Approximation 
,9]: asc = j:V nDk. This offers a good estimate of the 
scattering cross section. We have plotted the scatter- 
ing cross section over 175c- Plotting this quantity re- 
moves the energy dependence of the scattering. In Fig. 
[5] we have plotted both the total and inelastic (m chang- 
ing) cross sections over asc ■ The different energies are 
fe=10 (black solid line), 35 (red dashed), 53 (blue cir- 
cles), and 80 (green +). This estimate is given for the 
distinguishable case. In the case of bosons or fermions 
their a will oscillate about aEik @, EH- Notice that the 
elastic scattering never turns off even though there is no 
diagonal interaction at ol c /tt ~ 0.30. This shows that the 
scattering is made up from second-order processes; even 




FIG. 2: (Color Online) The total and inelastic a /asc as a 
function of a is plotted for four different values of Dk: 10 
(black solid line), 35 (red dashed), 53 (blue circles), and 80 
(green +). The inelastic data is when m is changed in a 
collisions, and in this case, it makes up a significant fraction 
of the scattering. 



though there is no diagonal interaction, there is a sig- 
nificant diagonal scattering contribution because of the 
off-diagonal channel couplings. The shape of this curve is 
that the total scattering rate dips and reduces to about 
40% of it original value then increases up to about 65%. 
The inelastic rate starts at zero and quickly climbs until 
a/ir > 0.2, but after that it only slightly increases. For 
the lower end of this regime (k = 10) scattering, there 
are still noticeable resonances in the scattering when a 
single partial wave makes up a large percentage of the 
scattering. 

With these simple estimates of the scattering magni- 
tude in hand, we now move to study the impact of the 
tilted polarization axis on the scattering for bosons or 
distinguishable dipoles, where there is the m = contri- 
bution relaying information about the short range scat- 
tering. For fermions, the short-range is strongly shielded, 
and only when one considers specific cases does the scat- 
tering become more involved. For this reason, we are 
more interested in the bosonic scattering and general 
long-range behavior and leave the case specific fermionic 
scattering for the future. 

As a first look at this scattering behavior, we look at 
the scattering rate as a function of po/D at three dif- 
ferent angles: a) a/n = 0.25, b) a/n = 0.35, and c) 
a/n = 0.50. In Fig. [3] a) the resonances which occur 
as po/D is decreased are narrow. Only a few exist be- 
cause there is a dipolar barrier to the scattering and thus 
the inter-channel couplings at short range must be strong 
enough to support the bound state. As the polarization 
angle is increased, the resonances become wider and more 
frequent. As the barrier is turned off and ultimately turn- 
ing into an attractive potential, we see more frequent and 
wide resonances. These resonances are much like the long 
range resonances seen in 3D [22T-|27j] . 

To study this energy dependence of the scattering fur- 
ther, we look at a particular angle and vary the energy. In 



4 




FIG. 3: (Color Online) The scattering rate, kao, as a function 
of po/D at three different angles: a) a/n — 0.25, b) a/n — 
0.35, and c) a/ir = 0.50 at Dk = I0" 2 (black). Notice that 
these resonances for small angles are narrow and fewer in 
number than the case where polarization is in-plane. In that 
case there are many resonances and they are wider in po/D. 



4 



ttf 







.-. ! 


. i 
- ji 
■ '. '[ 


■ i| 




II 

." . M 
■ ■ 1 t 


A Ai 






'. '.! 

: ■ 
. . 


" ; i 


■ V 

: j 


■'I 


■>. 


; - . ' 1 

■ i 1 

- .i 
-i 


■ f ( 




: : 
' '"j 


■ ■ 


.' i 
■ i 
; t 
. f 

; !■ 


; ; 
. i- 
■ r 

j ■ 


■ t 1 
i- 

I - j 
1 ■ 1 1 


1 

i 


j 

i- i 
' " / 1 


;■ 1 
i '■ } 
i ■ j 




j 
1 

' / 








y \ : 




' "■/ 

' / 
















T- 


- Dk=10" 3 
" Dk=10" 2 

- Dk=10"' 
■ Dk=l 


— 1 — 1 — 















FIG. 4: (Color Online) Energy dependence of the scattering 
rate, kao and the binding energies are shown as a function 
of po/D for a/ir = 0.35. The curves are for Dk = 10~ 3 
(green with circles), 10~ 2 (black), 10 _1 (dashed red) and 1 
(blue dotted). The binding energies are also shown for the 
same po/D. The cross sections and binding energy go to 
zero simultaneously and the strong energy dependence of the 
scattering resonance width. 
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FIG. 5: (Color Online) The scattering rate, kao, and a 2 Et, are 
shown as a function of a for many values of a. The different 
energies are: Dk = 10~ 3 (green), 10 -2 (black) 10 _1 (red) 
and 1 (blue). Cyan open circles are the effective range at 
each energy. The energy dependence of the scattering rate 
is seen simply as where the scattering rate peaks when ka = 
2e ' ~ 1.12. The binding energies times a 2 are also shown 
for all resonances, the blue dashed line is the universal limit, 
~ 0.63. Only for the Dk = 1 does the effective range not give 
a good estimate of the scattering rate. 



Fig. |4]we have plotted the scattering rate (a/n = 0.35) 
in the upper panel at four values of Dk: 10~ 3 (green 
with circles) , 10~ 2 (black), 1CP 1 (dashed red) and 1 (blue 
dotted). In the lower panel we have plotted the binding 
energy. There are a few points to be made here about 
the strong energy dependence of the scattering. First, 
the peak of the resonance shifts noticeably as the energy 
is lowered. Second, the width of the resonance becomes 
more narrow as the energy is decreased. Third, as the the 
binding energy goes to zero, the scattering rate goes to 
zero; this is in contrast to 3D. Fourth, for the \m\ > res- 
onances are very narrow, and they bind tightly as po/D 
is decreased. In these plots these resonance are hard to 
see because they are so narrow. They are most easily 
found by looking at the binding energy where there are 
steep lines. 

It is worth commenting on the relationship between 
the scattering length and the cross section and how these 
two quantities relate to the binding energy. First, for 
identical bosons in 3D, the cross section is a ~ 87ra 2 but 
for a 3> k it saturates at 8ir/k 2 . In the large a limit, the 

binding energy goes as ^— 2 ^ ^ , and the maximum of a 

corresponds to the maximum of a. 2D is very different. 
Consider when a — > oo, and the scattering cross section 
goes to zero. This is most easily seen as from the effective 
range expansion: cot(<5) oc ln(fca) — > oo. This leads to 
5^0 and a ~ when a is very large. The maximum of 
the scattering cross section occurs at 8 = 7r/2, and this 
happens when ak — 2e~ 7 ~ 1.12. 

A way to clearly demonstrate the energy dependent 
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behavior is shown in Fig. [5j where we have replotted 
the scattering rate and binding energy as a function of 
a/D, not po/D. We have replotted all the scattering 
data (i.e. many different values of a) as a function of 
a/D for many different values of Dk: 10~ 3 (black), 10~ 2 
(red) 10 _1 (green) and 1 (blue). The binding energies 
are now plotted as a 2 E[>, this makes it so that when the 
energies become universal they go to a constant value 
of ^-0.63. Replotting the data this way, allows us to 
observe clearly several points. First, the scattering rate 
is maximum when ka = 2e -7 ~ 1.12; this is clear from 
the four different energies plotted. It is also clear that 
the 2D system has strong energy dependence, and that 
the scattering rate goes to zero in the large a limit. 

From this plot we can understand why there is such 
strong energy dependence to the width of the resonances 
as a function of po/D. a/D is energy independent once 
one is in the threshold regime and depends only on po/D. 
We know that the scattering rate is zero when a — > oo 
and that it is maximum when ka ~ 1.12. Therefore as 
energy is decreased, the maximum and minimum of the 
scattering rate approach each other in a/D (Fig. [5]) or 
p /D (Fig @]). 

The effective range expansion gives a very good es- 
timate of the scattering rate and therefore the phase 
shift at low k. For k=l (blue squares), we are leaving 
the threshold regime, and the effective range description 
breaks down. 

Moving to the binding energies, we are see a 2 Et, (black 
circles) converges to the universal value (blue dashed line) 
for a/D > 10 and when it is strongly system dependent. 
We also see that at small a/D the binding energies de- 
viate from the universal trend and a 2 Ei, widely varies. 
In this figure we have only plotted the binding energies 
which were numerically found. Going beyond a/D=100 
is both computationally challenging and only returns the 
universal binding energy. 



IV. MOLECULES 

In this system the molecules have widely varying 
properties depending on the polarization angle. To 
study them more closely we pick six angles to explore: 
a/7r=0.25, 0.275, 0.3, 0.35, 0.4, and 0.5. For a/n smaller 
than 0.2, very little variation in the scattering, and there 
are no bound states for the values of po/D we consider. 
The first important point is that for a given angle the 
properties are robust. This is demonstrated by obtain- 
ing the molecular energies and wavefunctions for the first 
two m = resonances for each angle. Then we determine 
the values of po/D at each resonance which result in a set 
of chosen scattering lengths. We pick 40 different scat- 
tering lengths between 0.1 and 100D that are found at 
each of the first 2 resonances for each of the 6 angles. 

This idea of using a/D as the control parameter was 
used in 3D dipolar scattering to study three body recom- 
bination of dipoles [28[ . That work showed that using the 




Scattering length (a/D) 

FIG. 6: (Color Online) The binding energies for the first and 
second m — resonances at a/ir= 0.25, 0.275, 0.3, 0.35, 0.4, 
and 0.5. The universal 2D binding energy at large a/D is 
shown as a blue dashed line. 



scattering length to characterize the 2-body system in 
the calculations, even outside the a~^> D regime, worked 
well at revealing universal characteristics of three body 
recombination. 

In hg. [5] we plot the binding energies of the molecules 
for the first (black line) and second (red circle) reso- 
nance for all values of a considered. For a/ir = 0.25 
the molecules are most tightly bound at small a/D and 
at a / = 0.5 are most loosely bound, as expected, al- 
though there is roughly 2 orders of magnitude difference 
in the binding energies between these two extreme cases. 
For the tightly bound case the binding energies plateau 
as a/D is lowered. This energy corresponds to the size 
of the hard sphere and therefore a minimum size of the 
molecule. In this case there is a strong dipolar barrier 
and it is the inter-channel couplings that form the at- 
tractive short range region where the molecule is found. 
In contrast for a/n — 0.5 there is an attractive dipolar 
term for the m = case and the molecules are only mildly 
multi-channel objects. This will be explained below. 

As a/D is increased, for all of the polarization an- 
gles, the system becomes more loosely bound, and they 
are strongly system dependent. Only for relatively large 
a/D, say 10, do the binding energies truly resemble the 
universal value. We have plotted the universal binding 
energy as a blue dashed line. 

We now consider the size and shape of the molecules. 
First in Fig. [7] (a) we look at the expectation value of the 
molecular size: (i/j\p\i/j) as a function of scattering length. 
All values of a are plotted for both the first (black line) 
and second (red circle) resonance. For a/ir = 0.25 and 
small a/D the molecules are very small, (p) ~ 0.05D; this 
is roughly the size of the hard sphere when the first bound 
state is captured. In contrast a/it — 0.5 the molecule 
is about (p) ~ ID even for a/D ^0.1. In the large 
a/D limit we find that the size of the molecule goes to 
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FIG. 7: (Color Online) (a) The size of the molecules are 
shown as a function of a/D for various angles for both the 
first (black) and second (red circles) resonance. The univer- 
sal form of the size is shown as a dashed blue line, (b) The 
partial wave population is shown (first resonance only) , for all 
six angles as a function of a/D. The extreme angles are shown 
as filled circles (a/n = 0.25) and dashed lines (a/n = 0.5) for 
partial waves m = (black), |m| = 2 (red) and |m| = 4 (blue). 



7r 2 e 7 a/32 ~ 0.55a (blue dashed) which was obtained from 
4>b(K a p). Again we find that the size of the molecules 
depends on the polarization until a/D > 10. 

Now we look at the shape of the molecules shown in 
Fig- El (b). This is done by considering the partial wave 
population n rn = (4> m \4>m) as a function of scattering 
length for each angle. The extreme angles of a/ir = 0.25 
and a/ir = 0.5 are shown as filled circles and dashed lines 
for m — (black), |m| = 2 (red) and |m| = 4 (blue). 

For a/ it = 0.25 and small a/D the molecules are 
strongly aligned along the polarization axis behind the 
dipolar barrier. This is why they are so small and highly 
aligned. This is seen by the fact that the largest contri- 
bution is from |m| — 2 (red circles) and \m\ — 4 (blue cir- 
cles) is nearly 10% of the partial wave population at small 
scattering length. This is from the fact that the dipoles 
are behind the dipolar barrier and the inter-channel cou- 



plings are the origin of the molecule. The anisotropy at 
small a/D is still true for a/ir — 0.5, but the m=0 con- 
tribution is nearly 80%. This strong contrast is from be- 
cause of the attractive dipolar interaction for the m = 
channel only requires a slight amount of inter-channel 
coupling to form a bound state. For the 3D case, the 
/ = molecule is made up of about about 60% s-wave 
and the rest is essentially d-wave [281 ] . 

It is important to notice that at large a/D the par- 
tial wave populations go to very similar values. The 
m = contribution dominates and all other contribu- 
tions become small. To better understand the shape of 
the molecules we plot the radial weighted molecular den- 
sities. 

In Fig. [8]we have plotted the radial weighted molecular 
densities: p\if>(p)\ 2 = \ J2 m "m(p)e lm0 | 2 for a/n = a) 

0. 25, b) 0.35, c) 0.5 (top to bottom) for a/D of (i) 0.1, (ii) 

1, (iii) 10, and (iv) 100 (left to right). These densities are 
generated from the second resonance. The first resonance 
wavefunctions look the same expect there inner hard core 
is larger and over takes the inner oscillations. The scale 
changes for the plots on the left. 

This plot clearly shows the change in both shape and 
size of the molecules as both a/D and a/ir are changed. 
Now starting at a/D = 0.1 and ev/7r=0.25 (ai) we see 
that the molecule is very small and highly anisotropic. 
The density is along the polarization axis, x, and very 
tightly bound against the hard core. In fact its spatial 
extend does not really extend beyond the hard core in 
the y direction. 

Then as the angle is increased, for a fixed a/D the 
size of the molecule and anisotropy is softened. Observe 
the change is scale in both (bi) and (ci). In (ci), the 
molecule is larger (~ ID), and still aligned along the 
polarization axis. Additionally see the contrast in (ai) 
and (ci) between the extend of the density over the width 
in y of the hard core. 

Now for (a), (b), and (c) consider increasing a/D. The 
size of the molecules gets larger and more isotropic. In 
(iv) for all angles the system is isotropic, except for small 
region near the hard core, but the bulk of the radial 
weighted density is isotropically distributed at large p. 
This is why the size and shape of the molecules are uni- 
versal in the large a/D limit. 



V. CONCLUSIONS 

In this paper we have studied the scattering properties 
of the a pure 2D dipolar system when the polarization 
can tilt into the plane of motion. We have shown how 
the tilt angle impacts the scattering in both the thresh- 
old and semi-classical regimes. We then studied the char- 
acter of the scattering resonances generated by altering 
Pa/D or electric field as a function of the polarization 
angle. We found that at small angles the systems gained 
bound states which produced narrow resonances. This is 
because of the dipolar barrier. We also found when the 
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a)a/7i=0.25 




-50 50 

x/D 

FIG. 8: (Color Online) p\tp{p)\ 2 for: q/tt = a) 0.25, b) 0.35, c) 0.5 for a/D (i) 0.1, (ii) 1, (iii) 10 and (iv) 100. The contours 
indicated are drawn ever 20% of maximum value, with the largest 80% contour being drawn as red. The solid grey circle in 
the middle is the hard sphere short range interaction. The scale changes for every plot on the left (i, ii), in contrast the scale 
is the same for the plots on the right (iii, iv) when the system is in the large a/D regime. 



polarization is entirely in the plane of motion, the res- 
onances are frequent and wide because of the partially 
attractive potential. 

We studied the molecular system generated the tilt of 
the polarization. We showed that at large a/D the sys- 
tem have a universal shape independent of polarization 
angle, but at small a/D we have found that the molecules 
have a wide range of properties which strongly depend on 
polarization angle. Future work on this topic will be to 
more fully consider a realistic system, to include the ef- 
fect of confinement, fermionic dipoles and layer systems. 
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